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Abstract 

Recent two-photon photo-luminescence experiments give accurate data for the ground 
and first excited excitonic energies at different nanotube radii. In this paper we compare the 
analytic approximations proved in |CDR| . with a standard variational approach. We show an 
excellent agreement at sufficiently small radii. 

1 Introduction 

Recent experimental results on carbon nanotubes using two photon photo luminescence |Sc02j . 
[W] . |MPR| reveal the energy levels of the excitons, especially the ground and first excited states, 
and point out the dependence of these energies on the radius of the nanotube. As it is, this 
technique appears to be a promising way to sort out nanotubes. But on the other hand, theoretical 
results seem to require heavy ab initio calculation like in |MPR| [CBRMj . to cite the most recent, 
in order to find the absorption peaks due to excitons. Nevertheless, a simple formula for the 
optical response based on excitons levels |HKj can give a good approximation. It has been pointed 
out [BGEHj [OTj that the exciton binding energy in quantum wires depends on the width of the 
wires by a relatively simple relation. This property is also valid for nanotubes [PlJ. In the first 
part of the paper, we outline a rigorous justification for this latter fact and write a simple analytic 
formula for the energy levels of the exciton depending on the radius of the tube, based on the 
paper [CDRj . In the second part we compare our results with a variational numerical method and 
show a very good agreement between them. 

2 The exciton model 

As first suggested in [PI] . we deal with Wannier excitons (a rigorous justification of this procedure 
is in preparation [R]). We take as configuration space a cylinder of radius r and infinite length, 
space denoted by C r = M. x rS 1 , S 1 being the unit circle. The coordinates on the cylinder are 
(x, y) £ (R x rS 1 ) where x is the variable along the tube axis and y is the transverse coordinate. 

The two quasi-particles live in the Hilbert space L 2 (C r x C r ). We formally consider the Hamil- 
tonian 

H r = -h 2 ( ^ + + p± + p*) - V*{ X1 -x 2 , yi - y 2 ), (2.1) 
\Zmi Zm2 im\ Zm,2 J 

where 

V r (x,y) : = ^ = ' ei62 (2.2) 
eJx 2 + 4r 2 sin 2 (^) 
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(xi,yi) are the coordinates on the cylinder of the two charged particles, rrii their masses, and ej 
their charges. Here e is the electric permittivity of the material. The potential V r is the three 
dimensional Coulomb potential simply restricted to the cylinder. We justify the expression of V r 
by Pythagoras's theorem. The cylinder is embedded in M 3 . The distance p from one particle to 
the other in R 3 is: 

2 / v2 , a 2 • 2 fVi ~V2\ 

p ={x 1 -X2) + 4r sin I — — — I 

where |2rsin ( t,1 2 / /2 ) I is the length of the chord joining two points of coordinate y\ and y2 on the 
circle. 



2.1 "Separation of the center of mass" 

Due to the restrictions imposed by the cylindrical geometry, the usual separation of the center of 
mass does not work here. 

The standard separation with Jacobi coordinates only works for the longitudinal variable, and 
introducing M := mi + m,2 and p := miiri2 / (mi + TB2), we denote X := (m±Xi + ni2X2)/M and 
x := x\ —X2- For the transverse variable, we change to atomic coordinates Y = 7/2 and y = y\ — y2- 
Let us also define the effective Rydberg Ry* = pe 4 /2h 2 e 2 and Bohr radius a* B = h 2 e/pe 2 , where 
we set e = ei = B2- By a scaling, the new energy and radius will be expressed in multiple of these 
units. This gives us the Hamiltonian: 

H r = -^d 2 x ~—d 2 Y --d 2 ~-d 2 + —d y d Y -2V r (x,y). 

M 7712 M A* TO 2 

First, we can separate the partial center of mass with coordinate X. Second, since on the Y variable 
there are periodic boundary conditions, let us consider the orthonormal basis of eigenvectors of 



-d 2 



-dy= >. K\x r n )(x r n \ 

where 



00 

E 

n— — 00 



1 y n 2 
X r n (Y) = -j=e m ^ and E r n = ^, n G Z. 

One can see that for small radii, the separation between different transverse levels of energy is 
high. A recent theoretical study based on ab initio calculations |MPR) shows that the probability 
density of the exciton is constant along the circumference. So it is reasonable to assume that the 
radius is so small that the system stays in the ground state of — d Y , where 71 = and the density 
of probability is constant along the circumference. (Note that this has been rigorously proved in 
|CDR| ). As a consequence, we can approximate the eigenfunctions ip of H as: 

ip(x, y, Y) = <f>(x, y) ■ Xo(Y) = (j>(x, y) 



After this restriction to the lowest transverse mode, we only have to study the following operator 

iT = -±d 2 -l%-2V r (x,y), 

p p 



which is two dimensional and acts on L 2 (C r ). 



2.2 An effective one dimensional operator for the low lying spectrum 

It will turn out that the limit r — > is too singular and H r does not have a "nice" limit. It is 
suitable at this point to introduce the quadratic form associated to H r , defined on the Sobolev 
space ^(Cr): 

t^A) - -(d x ip,d x 4>) + -(d v i>,d y <p) - 2(^V r (x,y)i;, jV r {x,y)<l>). 



2 



Reasoning as in the previous subsection, a good approximation of the behavior of eigenfunctions 
along the transverse variable when the radius is small is given by the ground state of the free 
Laplacian with periodic boundary conditions. If we restrict the above quadratic form to functions 
of the type: 

</)(x, y) = tp(x) ■ — , (2.3) 



we have: 



'2nr ' 



tjj{(j>, <j)) = — — (d x (p, d x tp) - — I I ; = = ip(x)Tp(x)dydx 



Ivrn irr J_ x J_ nr / 2 + ^ ^ J 



(d x ip)(x)(d x tp)(x)dx — 2 / Vc S {x)(p(x)ip(x)dx 



oo 



where 



i r r 

AIU J— ITT 



-.dy. 



+ 4r 2 sin 2 if- 

2r 



For the sake of simplicity, we will put fx = 1. Let us introduce what will be our effective one- 
dimensional comparison operator: 

K S --=-^-W; s (x). (2.4) 

We now have reduced our problem of two particles on a cylinder to a one dimensional problem 
describing a particle interacting with an external potential. A complete and detailed mathematical 
justification of these steps can be found in [CDRj . 

2.3 The one dimensional Coulomb Hamiltonian 
2.3.1 Approximation of H^ s 

We now try to approximate as accurate as possible the eigenfunctions and eigenvalues of H^ s when 
the radius is small. Using the notation Y r {x) = ^/ x i +ir i i h was shown in |CDR| that if r < 1 and 
/, g are two smooth functions 

(f,V c r s g) = (f,Y r g) +ln(4)/(0)^(0) + O(r 4 / 9 )||/|| wl(K) || ff || wl(R) . (2.5) 

One can recognize in Y r a form of the regularized one dimensional Coulomb potential. In the 
following, our approach will be quite similar to the one of Loudon [L] in the sense that we let 
the regularizing parameter tend to zero. Note though that our parameter has a clear physical 
interpretation being given by the radius of the nanotube, and it is not just an artifact as in 
Loudon's case. Hence this will give us the exciton energies as functions of r, and allow us to 
estimate the errors we make using this approximation. 



Define the quadratic form Co(., .) bjo 

Coif, /):=-/ M*) ' (IfH'Wdx + I ln(-x) • (\f\ z )'{x)dx 



oo r 
2\l i 





° M-x) ■ (l/l 2 ) • (\f\ 2 Y(x)dx 

^ -L . \f( x )fdx. (2.6) 



1 Notice that this definition of Co differs from the one of [CDRJ in which there is an additional term: 21n2|/(0)| 2 = 
f_ ln2 ■ (\f\ 2 )'(x)dx + J °° ln2 ■ (\f\ 2 )'(x)dx. This term is here put together with the other term depending on 
the function at zero. 
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Note the last equality is obtained by an integration by parts, and holds for all e > 0; also note the 
appearance of a Coulomb potential in one dimension. Again from |CDR| we have 

(f,Y r g) = -21n(r)/(0)^0) + C (/,5) + O(r 4 / 9 )||/|| wl(R) || 5 || wl(R) . (2.7) 

Now let us define a new "potential" via a quadratic form 



(f,V C g) := -21n(-)/(0) 5 (0) + C (/, 5 ), 

which is close to when r is small and is exactly the Coulomb potential when we look away from 
the origin. We will see in the following that this particular potential gives a solvable eigenvalue 
problem for the associated Hamiltonian. 

2.3.2 Boundary conditions of the Coulomb Hamiltonian 

The operator He we will now consider is given by its associated quadratic form defined on ?< 1 (R): 



tcty,4>) ■= / <f>'{x)^{x)dx + 2 21n(-)0(O)V(O)-C o (^V) 

One can recognize the kinetic energy in the first term and the potential Vc in the second term. 
The general theory of closed symmetric quadratic forms gives us the existence of an associated 
operator He defined by 

t c {ct>,i)) = {cP,H c i>}, (2.8) 

whenever ip is in the domain of He- The difficulty here is that Vc is not a usual Schrodinger-type 
multiplication potential, but due to various Sobolcv embeddings it turns out that if i)j is in the 
domain of He then ip" is square integrable outside the origin, and we still have: 

(Hd>)(x) = - x^O. (2.9) 

Now if ip is an eigenfunction of He corresponding to an energy E, then it obeys the differential 
equation: 

-ij»( x )-WW=EiP(x), x^Q. (2.10) 
\x\ 

In order to get the behavior at the origin of the eigenfunctions of He, we integrate by parts 
and use (|2.10p . For e > we have: 



(j)'(x)ij'(x)dx (2.11) 



e)ip'{-e)- / <j>(x)^' {x)dx - <I>{e)^ {e) - / (j){x)ip" {x)dx 



+ J (/)' (x)ip' (x)dx 

where the last integral will converge to zero as e goes to zero. On the other hand, 
C (<p,yj) = [ ln(-x) ■ (d x {W)){x)dx - [ In(x) • (d x (^))(x)dx 



+ ln(e) (cj){e)-p{e) + cj)(-e)^{~e)] + / — ■ (f>{x)iP(x)dx. (2.12) 

V > JR\[- E ,e] Fl 

Then adding ([2~TTj) with (|2~T2|) . using (f2TT0|) and letting e tend to zero (see }CDRj for technical 
details) we have: 

1 21n(^(0)-Me)(2V(0)) 



lim t C (& VO = (0, #V>) + 2 lim 0(0) 

£— >0 s— *0 



2 v 2' 
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Now using (|2.8p and the eigenvalue equation Hctp = Eip we get the boundary condition: 



lim 



^(-e)-^(e) 



21n(-)V(0)-21n( £ )V(0) 



= 0. 



(2.13) 



2.4 Eigenvalues and eigenfunctions 

We now have to solve the equation 



(2.14) 



with the boundary condition (|2.13p . Similarly to Loudon in (Lj, we introduce a dimensionlcss 
quantity a and the change of variables 



1 , a 

E = - and x = —z, 

a z 2 



then we obtain 



z^O. 



(2.15) 



(2.16) 



The solutions are known for z > and z < 0, see |AS[ chap. 13] for example. The second thing is 
to see what condition at z = should the eigenfunctions -0 obey. If we scale (12.131) , 



lim 



+ a(hxr — ln(ae))tp(Q) 



0. 



(2.17) 



We only give here the results. Details of calculations can be found in |CDR| . If we take the only 
square integrable solution, we have two cases: If a = N is a positive integer, then the eigenstates 
are the odd functions tp na p with associated eigenvalues E notP , where n a = N + 1, 



^n aP (z)=e i^z— ==L] v _ 1 (|z|), 



2N 



E =-± 



where L 1 N _ l is an associated Laguerre polynomial. Notice that these energies are independent of 
the radius of the tube. If a is not an integer, the eigenstates are the even functions ipn a s with 
associated eigenvalues E riaS , where n a is the smallest integer larger than a and: 



C a W a 



i(|z|) = C Q |z| e -^l[/(l-a,2,|z|) 



1 



where C' a is a normalizing constant, W is the Whittaker function and U is the confluent hyper- 
geometric function or Kummer function of the second kind. We denote with T(z) and ^(z) = 
T'(z)/T(z) the usual gamma and digamma functions, and we get from (|2.17|) that for even solutions 
a must satisfy the relation: 



*(1 - a) + 27- 



In a + bar = 0. 

2a 



(2.18) 



From this relation, which contains an implicit expression for a(r), one can deduce several important 
facts. For all integers N and for a in between N and N + 1, there is only one solution of (|2.16|) 
satisfying the boundary condition. Furthermore, the energies associated with the non integer a 
tend to their closest lower integer when r tends to zero. A special case is the one of the ground 
state Ei s of the exciton which tends to minus infinity as the radius tends to zero. The behavior 
for small r is 

E u r =°-4(lnr) 2 . 
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•E 4p =-l/9 
E 3s (r) 
■E 3p =-l/4 



'E 2p =-1 



-EJr)~-4(ln(r)j 

Figure 1: Energy of the bound states of the exciton. Energies are expressed in multiples of the 
effective Rydberg of the exciton. 

Notice that equation (|2.18|) is exactly what one gets by explicitly calculating the condition 
requiring the derivative of even states to vanish at the origin, relation (3.22) in Loudon's paper, 
and replacing his parameter a by r • ao/2, for small r. We can also compare with the result 
from [BGEH] where Banyai et al. found numerically a = Q.'Sao ■ r for the exciton problem in a 
quantum wire of radius r. The energies associated to the exciton are drawn on figure [T] and we 
calculated numerically the lowest four eigenvalues with respect to the radius on figure [3J along 
with a comparison of the ground state given by the variational method of |Plj . On figure [4] a zoom 
was made around the first and second excited states. 

One can notice that the odd states possess a constant energy, independent of the radius. This 
is due to the fact that the odd states vanish at zero where there is the singularity of the potential. 
Indeed, the potential becomes, for an odd function /: 

(f,v: e f)~ Co (/,/), 

which is independent of r. 

3 Variational approach and comparison 

The variational method operates by minimizing the energy of a trial function and is therefore 
usually applied to approximate the ground state. However, by restricting the trial function to 
forms that are orthogonal to the ground state, low lying excited states can be obtained variationally 
as well. In Ref. [PI] , the variational method was applied to the ground state and in Ref. |P3| 
a similar approach was applied to calculate the 2p oscillator strength of interest for two-photon 
absorption. The 2p state is especially important because this state is used for excitation in two- 
photon fluorescence measurements |Wj . By recording the energy of photons emitted from the 
lowest (Is) exciton, a direct measure of the 2p — Is energy difference is obtained. In turn, the Is 
exciton binding energy (i.e the Is excitation energy measured relative to the band gap) can be 
derived if a reliable model for the exciton energy spectrum is invoked. 

In the present work, we wish to compare results of the relatively complicated variational 
approach to the straight-forward and analytical Coulomb model. Hence, in the following we 
briefly explain the reasoning behind our variational estimate of the lowest excited (2p) state. In 
practical applications, the trial function must be sufficiently simple that calculation of the energy 
is manageable. This implies that relatively few adjustable parameters should be considered. At 
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-0.4 




Nanotube radius r [a* ] 

Figure 2: Variational 2p state energy as a function of nanotube radius. Inset: the r-dependence 
of decay lengths k and q of the trial wave function. 

the same time, a certain flexibility is required to provide reasonable accuracy. A useful strategy 
consists in constructing trial functions so that they correctly accommodate the known solutions 
in limiting cases of the general problem. Thus, we are guided by the analytical solution for the 
plane, i.e. for nanotube radii much larger than the effective Bohr radius. In our units, this state 
is simply ip2 P oc x exp{— 2/3(x 2 + y 2 ) 1 ^ 2 }. On the other hand, in the extreme ID limit we expect 
ifi2p oc a;exp{— \x\}. To accommodate both limits, we consequently suggest the (un-normalized) 
trial form 

V2 P (x, y) = xexp{-(x 2 /k 2 + y 2 /q 2 ) 1/2 }, (3.1) 

where q and k are variational parameters to be determined by minimizing the expectation value 
of the energy . This expectation value is found as E2 P = (K — V)/N, where K, V and N are 
the kinetic energy, potential energy and normalization constant, respectively. The integrations 
are quite cumbersome and only K and N can be obtained in closed form in terms of Struve 
and modified Bessel functions. The potential energy V is evaluated numerically using Gaussian 
quadrature. 

The minimized energy as a function of radius r is illustrated in Fig. [2] The limiting values 
are i?2 P = — 1 and Ei v = —4/9 for small and large r, respectively, in agreement with the analytic 
solutions in these limits. In between these limits, the curve interpolates smoothly between the 
limiting values and the dominant correction at small r is —8(1 + 7 + ln(r))r 2 . 

In order to judge the usefulness of the different approaches used in the present work it is 
essential to determine the appropriate nanotube radius r in excitonic units, i.e. in units of a* B . 
An important point in this respect is that, in fact, a* B varies between different nanotubes. The 
relation needed for conversion is a* B = 0.529Ae//x, where e is the dielectric constant screening the 
Coulomb interaction and fx is the reduced effective mass. Whereas e may be assumed the same for 
all nanotubes, \i must be derived from the curvature of the band structure and, hence, depends on 
both radius and chirality of the nanotube. However, detailed studies [P2j show that a* B is roughly 
proportional to the nanotube diameter and as a consequence r (in units of a* B ) is nearly constant 
and given by r w 0.1a* B if e ~ 3.5. It should be noted, though, that in media with little screening 
a larger r/a* B is expected. The smallness of r means that approximations based on expansion 
around r — arc expected to be accurate. 

In figure [3] and HI a comparison of variational energies and the results of the Coulomb model 
is given. In both the Is and 2p cases, reasonable agreement between the two approaches is found 
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First p state £" 



2s 



2p 

Ground state E 
— Variational E, 
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Figure 3: Energy of the four lowest bound states of the exciton with respect to the radius. Even 
states were calculated numerically using (|2.18|) . Notice that the energy of the second s state is 
equal to the one of the first p state at zero. In red: graph of the ground state computed with the 
variational method on the cylinder from |P1| . 



around r ~ 0.1a* B . Also, in both cases, the variational result lies slightly higher than the Coulomb 
model. The ground state has deviated by more than 100% from the plane value for radii around 
O.la^. An error of less than 5% is seen between the curves at this point. Note that if the curve 
from the Coulomb model is lower in energy than the variational one, this does not means that 
the approximation is better. Although in the variational approach the lower is the better since 
the exact solution is always below the variational result, in the Coulomb model the exact value is 
somewhere around the approximation, bounded by an error bound. In this work the error bound 
was not calculated because, unfortunately, the compromises made to get a simple formula implied 
a difficult calculus to optimize on the bound, even by numerical computations. For the first excited 
state, the Coulomb model approximation for the energy is independent of the radius of the tube. 
This must be a good approximation for very small radii. However, this approximation becomes 
increasingly inaccurate as the radius increases, as the exact value should tend to the energy of the 
problem on the plane with energy equal to — 4/9a^ as do the variational curve. 



4 Conclusion 

The Coulomb model applied to excitons in carbon nanotubes demonstrates that the energy as- 
sociated with each even state decreases to the energy of its closest odd state as the radius tends 
to zero. Since each odd state energy is independent of the size of the tube in this approach, all 
energies associated to excited states stay in the range (—1, 0) (in effective Rydberg energy units). 
So only the ground state diverges when the radius gets smaller. This is confirmed by the behavior 
of the first variational excited state which converges to — 1 . The simple exciton model proposed 
in this paper and studied both by means of rigorous perturbation theory and by a variational ap- 
proach is a good starting point for attacking the difficult problem of electron-electron interactions 
in low dimensional structures such as carbon nanotubes. 
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Figure 4: Zoom around the energy of the 2p state showing both variational and Coulomb model 
results. In addition, the 2s state obtained with Coulomb model is included for comparison. 
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